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Abstract 

We calculate the thermodynamic properties of a collection of N small Josephson junctions cou- 
pled to a single-mode resonant electromagnetic cavity, at finite temperature T, using several ap- 
proaches. In the first approach, we include all the quantum-mechanical levels of the junction, but 
treat the junction-cavity interaction using a mean-field approximation developed previously for 
T = 0. In the other approaches, the junctions are treated including only the two lowest energy 
levels per junction, but with two different Hamiltonians. The first of these maps onto the Dicke 
model of quantum optics. The second is a modified Dicke model which contains an additional 
XY-like coupling between the junctions. The modified Dicke model can be treated using a mean- 
field theory, which in the limit of zero XY coupling gives the solution of the Dicke model in the 
thermodynamic limit using Glauber coherent states to represent the cavity. In all cases, for an 
iV-independent junction-cavity coupling, there is a critical junction number N above which there 
is a continuous transition from incoherence to coherence with decreasing T. If the coupling scales 
with N so as to give a well-behaved thermodynamic limit, there is a critical minimum coupling 
strength for the onset of coherence. In all three models, the cavity photon occupation numbers 
have a non-Bose distribution when the system is coherent. 

PACS numbers: 74.50. +r, 74.25.Nf, 85.25.Cp, 64.60.Cn 
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I. INTRODUCTION 



When a two-dimensional array of Josephson junctions is driven by an applied current, it 
can radiate coherently Experiments showing this behavior have emphasized current-driven 
arrays of overdamped Josephson junctions The radiated coherent power from such ar- 

rays has been predicted to be proportional to the square of the number N of Josephson junc- 
tions in the array s|. More recently, coherent emission from underdamped one-dimensional 
Josephson arrays coupled to a single-mode electromagnetic cavity has been experimentally 
studied . In this work, it was shown that no coherent radiation is emitted below a 

threshold number N c of junctions, but above this threshold the array can radiate coherently, 
with emitted power again proportional to N 2 . Such behavior had already been predicted 
much earlier, on the basis of an analogy between a one-dimensional voltage biased series ar- 
ray and a collection of two-level atoms coupled to an electromagnetic cavity [8( ; the analogy 
suggests that this radiation is the Josephson analog of the population inversion that leads 
to coherent emission in a laser. 

A simple model Hamiltonian to describe this coherent radiation, taking into account the 
quantum-mechanical nature of both the junctions and the cavity, was suggested recently 
In this paper, the ground state of the model Hamiltonian is obtained within a mean-field 
theory (MFT). In agreement with experiment, the MFT predicts that there is a critical 
threshold number N c of junctions for the onset of coherence at fixed coupling strength. 

The mechanism for coherent radiation from a Josephson junction array resembles that of 
superradiance in a system of N two- level atoms coupled to a electromagnetic field [h]]. The 



which describes the system of identical 



latter system can be treated the Dicke model 
two-level atoms in a single-mode radiation field. Emission and absorption within the Dicke 
model have been extensively studied 
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131 ] . The predicted response of this two-level 



atom/radiation system agrees qualitatively with that of an array of Josephson junctions 
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151 ] . It has also been shown JJjJ, that a modified Dicke Hamiltonian, which contains 



an additional term resembling a dipole-dipole interaction between the junctions, is a better 
approximation to the cavity-j unction system than is the original Dicke model. 

In contrast to the Josephson/cavity system, the Dicke model can be solved in the thermo- 
dynamic limit, i.e. iV — > oo, volume V — > oo, and N/V — > const [l6, [l7, [l8]. The solution 
yields a continuous transition from a normal to a superradiant state at a critical coupling 
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strength, for fixed temperature T = 0. 

In this paper, we extend the model 9] to finite T. We find, again within MFT, that there 
is a critical threshold number N C (T) for coherence at a sufficiently low T. If T is increased 
at fixed N > N c (0), there is a continuous transition from coherence to incoherence at an 
iV-dependent temperature T C (N). To test the MFT, we compare its predictions with those 
of the Dicke model [3] and of the modified Dicke model [l^. Our Hamiltonian does 
not map exactly onto these models, because the individual Josephson junctions have more 
than two quantum levels, whereas the Dicke and modified Dicke models assume two-level 
systems interacting with a single-mode cavity. Nonetheless, when the parameters of our 
model systems are such that the lowest two levels of the junction are well-separated from 
the higher levels, the MFT agrees well with the two-level model predictions. We also show 
that, when applied to the Dicke model, the MFT is equivalent to a coherent state expansion, 
an approach known to give the solution of the Dicke model at large N. 

The remainder of this paper is organized as follows. In Section [Til we describe the MFT 
for N Josephson junctions interacting with a cavity at finite T. In Section IIIIl we review 
the coherent state treatment of the Dicke model, present a MFT for both the Dicke model 
and the modified Dicke model, and show that, when applied to the Dicke model, the MFT 
is equivalent to the coherent state approach. In Section IIVI we give numerical results for 
all three models. In Section [V] we show how the Josephson- cavity model of Section [Til can 
be mapped onto the Dicke model when the Josephson coupling is small compared to the 
charging energy; we also discuss how the parameters of all three models must scale in the 
thermodynamic limit. Section IVTl presents a concluding discussion. 

II. MANY JOSEPHSON JUNCTIONS INTERACTING WITH A SINGLE-MODE 
CAVITY 

A. Model Hamiltonian 

An array of N underdamped voltage-biased Josephson junctions in a lossless electromag- 
netic cavity having a single electromagnetic mode of frequency u may be described by the 
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following idealized model Hamiltonian: 

N 

H = H photon + ^2 Hjj. (1) 
i=i 

Here 

H pho ton = M° ta + -) (2) 

is the photon Hamiltonian, a' and a being the photon creation and annihilation operators 
for photons having angular frequency u, which satisfy the usual commutation relations, 
[a, a)} = l,[a, a] = [a^a*] = 0. The Hamiltonian of the ] th Josephson junction can be 
expressed as 

Hjj = ]f (rij - n 3 f - J cos 7r (3) 

Here U = 4e 2 /C is the capacitive energy of the junction, e is the electronic charge, C 
is the junction capacitance, rij is an operator representing the difference in the number 
of Cooper pairs on the two superconducting islands forming the junction, rij is related to 
the gate voltage across the ] th junction, J = hl c /(2e) is the Josephson coupling energy of 
the junction, I c is the junction critical current, and finally jj is the gauge-invariant phase 
difference across the junction. 
Explicitly, jj may be written 

= f A dl, (4) 



% J j 

where <ftj is the phase difference across the junction in a particular gauge, and A is the vector 
potential due to the cavity mode (given explicitly below) in the same gauge, $o = he/ (2e) is 
the flux quantum, and the line integration is carried out across the junction. The operators 
rij and <f>j are canonically conjugate and satisfy the commutation relation [rik,4>i\ = idki, 
which is satisfied if we use the representation rik = i-^- 
In the Coulomb gauge, V • A = 0, A can be expressed as 



y^(a + at)E(x), (5) 



where V is the cavity volume, and E(x) is proportional to the local electric field of the cavity 
mode, normalized so that f v |E(x)| 2 rf 3 a; = 1. Introducing a coupling parameter 



we may rewrite 7 3 - as 

lj = <t>j -9j{a + <J)- (7) 

Eq. (jSJ) suggests that typically gj oc 1/ for a given mode (provided that the cavity s/iape 
does not change as the volume increases. 



B. Mean-field approximation 

We now develop a suitable mean- field approximation for the Hamiltonian (jTJ), for both 
zero and finite T. We consider only the case of identical Josephson junctions, so that all 



gj = g and rij = n; the extension to non-identical junctions is straightforward 19] . We also 



assume that the coupling parameters g are weak. In this case, we can expand the cosine in eq. 



([3]), retaining only the term of first order in g(a + at), so that coS7 3 - ~ cos <pj+g(a + at) sin 
Within this approximation, the only part of the Hamiltonian that depends on both cavity 
and junction variables is 

Hint ~ ~gJ {a + at) ^ sin fa. (8) 

i 

The eigenvalues and eigenfunctions of H can now be found if we make the following 
mean- field approximation for Hi nt : 

Hint ~H™t= ~ gJ{a + at) ^(wifa) 

3 

— gJ(a + at) sin fa 

3 

+ gj(a + a ] ) ^ ( sin 4>j)- ( 9 ) 

3 

Here (...) denotes a canonical average at temperature T with respect to the mean- field 
Hamiltonian H m . H m is now given by 

H m = H™ hoton + B?j + H ?i (io) 

3 

where 

K hot(m = + ~) - gj(a + at) ^(sm^), (11) 

3 

HJ] = — (rij — n) 2 — Jcos0,- — gJ{a + at) sin fa, (12) 
J 2 
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and 

H™ =^J(a + a t }^(sin0 J ). (13) 

j 

Evidently, the first term depends only on the photon variables, the second is a sum of 
single-junction terms, and the third is simply a c-number. 

We now introduce the variable Xj = (sin </>,,•) = A, since Xj is independent of j. In terms 
of A, we may rewrite the photon term (fTTI) as 

Khoton = M« f « + \)~ 9J(a + a^)NX. (14) 

This is simply the Hamiltonian of a displaced harmonic oscillator. Its eigenvalues E n 
and normalized eigenfunctions ip n (x) are just E n = hw(n + |) — J 9 Jj A and ip n {x) = 
( x \ e -ip(x)/n\ n ) = e -^-{*)) 2 H n (y/^ix - (x))). Here H n is a Hermite polyno- 

mial, p = iyp^-{cft — a) is a momentum operator, (x) = 9J ^ X \[^t~ l * s ^ ne m ean displacement, 
and we have used the relation (a + a)) = 2J t ■ The canonical partition function corre- 
sponding to H™ hoton is just Z™ hoton = 2 ^^/ 2 , where £ = #JNA and /3 = l/{k B T). 

Next, we consider the Josephson junction Hamiltonian Hjj [eq. ( 1X21) ] . With the defini- 
tions K(A) = J^l + A^^j^y and = e -^-") M (0 - a), where tana = ^0^, 
the Schrodinger equation for the junctions, Hjjip((f>j) = Eif)(<f>j), reduces to the standard 



Mathieu equation 



20J: 

^ 2 + (a - 2gcos2^)yfa) = 0, (15) 

where Vj = (<f>j — a)/2, y(vj) = u((</>j — a)/2), q = —4K(X)/U, the characteristic value of the 
Mathieu equation is a = 8E/U, and we have used the representation rij = i-^- The allowed 
eigenvalues are determined by the condition that ip(<f)j + 2ir) = xjj(<pj), or equivalently, 
that y(vj + 7r) = exp(2m7r)y(t> j). The allowed solutions y v {vj) are therefore the Floquet 
(Bloch) functions of Vj, with Floquet exponent v = 2n + 2k, where k — 0, ±1, ±2, .... The 
corresponding eigenvalues of Hjj are labeled by the quantum number v = 2n + 2k and the 
parameter q, and may be denoted E(y = 2n + 2A;; g). For < n < 0.5, the lowest eigenvalue 
corresponds to k = 0, followed in order by = —1, 1, —2, 2, ... Including only these Floquet 
solutions, we can formally express the junction partition function as 

Z™. = e - pEi2n+2k > q) = Zf, (16) 

fc=0,±l,±2,... 

where the last identity holds for identical junctions. 
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We now determine the properties of the junction-cavity system within MFT. At T = 0, 
the system is in its ground state, and the approximate ground state properties can be 
obtained analytically as shown in Ref. For —0.5 < n < 0.5, the ground state energy is 

„ frw igJNX) 2 U . n . 

E 9 (X) = -y + 1 ^ ^ + ^-a(i/ = 2^; g), (17) 

where a(z/; g) is the eigenvalue of eq. (fl5l) corresponding to characteristic exponent z/ and 
parameter q. If we use the approximate analytical expression [9] 



a(f = 2n; q) ~ 4 I 1 - 4n 2 - */ (1 - 4n 2 ) 2 + — ) + An 2 , (18) 



we obtain 



, > Hoo (gJNX) 2 U , 



\ 



AJ 2 

1 + w 



1 + 4 



g^JNxY 
h/jj J 



Un 2 

+ (19) 



where U = U(l — An 2 ). A is determined by the condition = 0, which leads to 

hio \ 2 (. u 2 



^ t = °^ 1 -{W7n) [ l + lr-)- (20) 

Because the right hand side of this equation must be non- negative, the critical junction 
number N c (0) for a non-zero A at T = is 

N c (T = 0)~J^^U 2 + AJ 2 . (21) 

When iV < N c (0), X = corresponds to a minimum of the energy, but when N > N c (0), 
X = is a local maximum; the energy minimum occurs at A ^ 0. 

More generally, the exact solution for N C (T = 0) can be calculated from eq. (TlTI) . supple- 
mented by the condition = 0. The result is 



. J(a'g 2 JN) 2 - (huj) 2 , . 

where a'(A) = da ( u -^ n ' q '> 4y(A) . The corresponding critical number is 



where a'(0) = a'(A)| 



A=0- 
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For T / 0, because H m is the sum of several commuting terms, the total partition 

■m I yva 
photon \ J 



function Z m = Z™ hoton (Z™) N Z™ , where Z c = exp(—@H™). The corresponding Helmholtz 



free energy F m is 

F m = -kT In Z m = F™ hoton + NFJ 1 + F™, (24) 

where F™ oton = fcTln (2 sinh (/^/2)) - (gJ ^ A)2 , Ff = -k B T\nZJ, and F c m = 2 (gJ ^ A)2 . 
When the coherence order parameter A ^ 0, the Helmholtz free energy, for fixed g, is 
quadratic in the number of the junctions N. This quadratic dependence is a hallmark of 
the coherent state. 

The actual value of A is obtained from the Helmholtz free energy, using the condition 

-2* -a 

We have obtained A by solving eqs. (!24l) and ( J25l) self-consistently. These equations may 
allow for several possible values of A, of which we choose that value which gives the lowest 
F m . 

In all the above discussion, we have assumed implicitly that g is independent of N. The 
expected behavior when g depends on N is discussed below. 

III. DICKE MODEL AND GENERALIZED DICKE MODEL 
A. Model Hamiltonians 

In the previous section, we described a simple mean-field approximation for the statisti- 
cal mechanics of the junction-cavity system. This approximation includes all the junction 
levels, but treats the junction- cavity interaction only approximately. We now describe an 
alternative approach, which retains only the two lowest energy levels of each junction. In 
this case, the Hamiltonian reduces to the well-known Dicke model of quantum optics. The 
Dicke Hamiltonian ll[ is a simple model describing the interaction of N two-level systems 
with a single harmonic oscillator mode. It can be written (omitting the cavity zero-point 
energy) 

H mck * = hwa^a + + e*av; + . (26) 

3 

Here ej > is the energy level splitting of the j th two-level system at the ] th junction, and £j 
is a parameter characterizing the strength of the coupling between the harmonic oscillator 



and the ] th two-level system. The quantities a) and a are raising and lowering operators, as 
above. The quantities a J a are Pauli spin-1/2 spin operators and satisfy [u^Co] = lib^o, i, 
where a, (3, and 7 are cyclic commutations of (x,y,z). In order for the thermodynamic limit 
to exist, we must assume that £j oc 1/yN for large N. 

Besides the Dicke model, we a 
sion of effective two-qubit model 



so consider a modified Dicke model [15], which is an exten- 
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231 ] to the case of N coupled two-level systems, 

H MDicke = huja \ a + CL^i + Z j( JaL + $aai) + £ Sl jk {o{a k _ + ota*), (27) 

i ^ ' 0'*> 

where the last sum runs over all distinct pairs jk (i. e., not including j = k). The last term 
in eq. (l27l) is an effective direct junction-junction interaction. As discussed in Ref. [15], the 
Hamiltonian (1271) generally gives levels in closer agreement with the Hamiltonian (TjQ) than 
does the pure Dicke Hamiltonian f l26l) . when there is more than one photon excited in the 

n 

cavity. A simple derivation of this term is given in Ref. [15] . 

In order for the thermodynamic limit to exist in the modified Dicke model, we require 
not only that ^ oc 1/yN, and but also that Vtj^ oc l/(N — 1), as further explained below. 



B. Statistical mechanics of Dicke model using Glauber coherent state expansion 

The thermodynamics of the Hamiltonian (f26|) can be calculated in the limit iV — > 00, 
using a product basis consisting of the Glauber coherent states \a) for the photons and 
eigenstates of a l z for the two-level systems [nj, Q]- In this section, we briefly review this 
solution. 

The states \a) are eigenstates of the lowering operator, i. e., 

a\a) = a\a). (28) 

where the eigenvalue a is generally complex, since a is a non-Hermitian operator. The 
eigenfunctions satisfy the completeness relation (1/tt) dRe(a)dIm(a)\a) (a\ = 1, 

the integral running over the entire complex a plane. 
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In terms of this basis, the partition function of the Dicke Hamiltonian takes the form 

-? HDlcke (29) 

\a)\oi...a N ) 



J Dicke 



Tre 



E •• E 

o-l =±1 a N =±l J -°° J -°° 



oo rco 



dRe(a)dIm(a) _ 3H Di a , 
(<Jx...a N \(a\e 



OO /"OO 



dRe(a)dIm(a) 



7T 



-/3fiu>[a| : 



oo J — OO 



where 



7T 



n e< 



>o-,=±l 



1 



(30) 



(31) 



In order to evaluate the sums, following Refs. [l7] and 18|, we expand the operator 
exp(—(3H Dlcke ) in a Taylor series, assume that a/y/N and a) /y/~N exist in the limit N — > oo, 



and finally, within individual terms of the Taylor expansion, interchange the order of the 
double limits as follows: 

R / 7-7- T)irlce>\<r R 

lim lim 

N^oo R—>oo ■ 



-f3H 



r=0 



lim 7 lim 

R-*oo ^— ' JV->oo 
r=0 



Dicke\r 



V ! 



(32) 



Each sum in eq. flHDl can easily be evaluated, since it is just the trace of the operator 
exp(—(3hj), and the resulting partition function can be expressed as 



J Dicke 



dRe(a)dIm(a) 



7T 



2 / rdre-^ r2 I 2 cosh ^ « / 1 + 



(33) 



where we have introduced r = lad and written f°° f 

1 1 j —00 j - 



v r°o dRe(a)dIm(a) 

OO J —OO 7T 



2 J °° rrfr. 



To complete the evaluation of the free energy, we make the changes of variables = 



and y — jr. Then Z Dicke. takes the form 



Z Dicke = N I dyexp [N<j)(y)] 
'0 



where 



{-fihuy + 1 J> (2 cosh 



(34) 



(35) 



The last integral can be evaluated accurately using Laplace's method [24J. This method 
makes use of the fact that, if N ^> 1, the integral should be dominated by values of y near 
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the maximum of <p(y), which is determined by the condition 



<j>'{y) = 0. (36) 
In our case, 4>'(y) = —j3huj + i £\ — ^ 3 13 = tanh 1 + so ec L- (|36|) becomes 



- = -r^tanhf^Y (37) 



where rjj = y 1 + — and lies in the range 1 < rjj < oo. For large (3, tanh (^rVj^J ~^ 1- 
Thus, when -r- > j, ,■ — > <K?/) is maximum at y = 0. On the other hand, when -j- < 

£'2 

t~ ■> the allowed solutions to eq. (136]) depend on the value of (3. When (3 is smaller than 
a critical value (3 C given by 

then again 0(y) is maximum for y = 0. However, if /3 > /3 C , there is a non-zero solution ?/o 
for y determined by the equation 



Therefore, we can discuss the statistical mechanics of this model in three different regimes: 
(i) E; S ^ X5 (**) E, lr > X and T > T - and (***) E; if: > x and T < T c , where 

J 3 J 3 J 3 

T c is determined by eq. (|38|) . In the regimes (z) and (zz), the free energy is given simply by 
Jim ^ = -k B T4>{y) y=Q = -k B T± £ In (2 cosh ^) (40) 

3 

and the moments of the photon occupation number by 



t \ k 

> = (y k ) y=0 = 5 h0 , (41) 



a'a 
~N 



where the last result is obtained once again by using Laplace's method. On the other hand, 
in regime (Hi), 

F 

lim T7 = - k BT(p(y) y=yo (42) 

N— >oo iV 



11 



and 

({^-) k ) = (y k )y=y ^y k o (43) 
where yo is determined by eq. (1391) . 

Assuming that all junctions are identical, 6j —>■ e, £j — > £, Xj — > A, the conditions for the 
critical junction number and the critical temperature become 

nrDicke _ ^ 6 (AA) 
fcflT «cfce = * L (45) 

& 2tanh- 1 (-^ 
Furthermore, the moments ((a^a) fc ) = iVyo are obtained from 



hue I 16iV^ 2 y , /9e A 16iV£ 2 ?/o 



4^v 1 + ^ = tanh TV 1 + ^' (46) 

C. Statistical mechanics of modified Dicke model using a mean-field approxima- 
tion 

Next, we consider the modified Dicke model, eq. (1271) . The coherence transition can be 
obtained if we make the following mean-field approximation: 

aV ~ (aVl + ^iy-^Xat), (47) 



aa 



j 



o\o k _ ~ 



(a)4 + (4)a-(4)(a), (48) 
(4)^ + (^)4-(4>(^>, (49) 



olo\ ~ (oL)o% + (o*)ot - (aL){a k + ). (50) 

With the additional assumption that the £j's are real, the Hamiltonian (|27jl separates into 
a sum of three terms as follows: 

Tj M Dicke tjM Dicke , \ ^ tjAI Dicke , tjM Dicke 

H = h photon + 2^ H Jj + h c ) I 51 ) 

i 

where 

#S= e = ^(ata+^) + ^0((^)(« + « t ) +<<^>(a - a f )), (52) 

i 

= -E^« at x^> + < a ><4» -E^+x^)' ( 54 ) 
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and we have used the relations (a J + ) = (a J x ) + i(cr J y ) and {at) = (a J x ) — i{cr J y )- 

The free energy associated with each term in the above Hamiltonian can be evaluated 
separately. H^ t l ^ e is the Hamiltonian of a harmonic oscillator displaced in both momentum 
and position space. Introducing the operators 

« = y^(a + a t ), p = i\J^-{a ] - a), (55) 

we find that H^t<^ e can be rewritten as 

1 . N o moj 2 



tt M Dicke _ 1 / „ \2 / n2 

^photon - ^ yp - Po) + ~y~ y u - u o) 



(56) 



where po = sJ^Ylji^ylij an d u o — —\J j^s Ylj^D^j- The corresponding eigenvalues are 

^ ( „ ) = M „ + ^)- Sa^>g>±<^ . (57) 

Also, from the fact that (p) = p and (u) = u , we obtain (a) = — ^' 3 ^~\ ip)) = 



hi) 

jjMDicke j g foe Hamiltonian of a collection of non-interacting spin-1/2 particles in an 
applied effective magnetic field (which is not parallel to the z axis). The two eigenvalues of 

jjMDicke are reac [ily f ounc l to be 



jpMDicke _i_ 



e '- + 4(0(at)+^^K)) U^ + ^^iot)). (58) 



2 7 V ^2 

fc^j / \ k^j 



Finally, H^ Dtcke is just a c-number whose expectation value is just E^ Dtcke = H^ Dlche . 

For N identical junctions, €j — > e, £j — > £, and — > fi. Then, after some algebra, one 
finds that the ground state energy can be written in terms of a single expectation value 
K> = [K> 2 + (^) 2 ] 1/2 - The result is 



(59) 

(crj_) is again determined by the requirement that E^ rD ' tcke ((a±)) be a minimum with respect 
to (cr±), which leads to 



1 / hue 



<^> 2 = 1 - i Up-^-dm I • (60) 

13 



The condition (cr±) 2 > leads to the critical number of junctions 

h 1 2n 

N MDicke = ± e_ (Q1) 

^? J- 2€ 2 

above which (<t_l) 2 is non-negative. When £1 = 0, this critical number exactly corresponds 
to the critical number obtained in Section IHIi Thus, for the Dicke model, this MFT yields 
the same critical junction number as obtained from the coherent state analysis. 

At finite T, the properties of the modified Dicke model are obtained from the Helmholtz 
free energy p MD%cke . An analysis similar to that at T = again allows F MDlcke to be written 
as the sum of three terms, which for N identical junctions may be written 

F MDlcke = ~k B T In - Nk B T ln(Zf Di ^) + E%™%, (62) 



and 



where Z$™? = j-^; Zf D ^ = 2 cosh + 4<<r ± >* (^fi - "gfj 

^constant = ~ " 'h^^ ~ ^m^j ' wn i cn also includes constant contributions from photon 
and junction terms. 

As at T = 0, the optimal value of (a±) at finite T is obtained by minimizing p MD%cke 
with respect to (cr±), which leads to the following relation for (er_i_): 



Humeri 

tanh 



7?e ' 



(63) 



where rj — *J 1 + 16i J 2( -^J)2 L ^ ^1 — ^i^p"^ ) ■ ^ ne CI "itical temperature T^ IDlcke for this 
modified Dicke model is again determined by the requirement that (cr±) 2 > 0, and is given 
by 

kgT M Dicke = 1 ^_ (64) 



2 tanh 



1 / hjjt 



4JV£ 2 (x hu(N-ljQ \ 

V 2JV^2 y 



k T M Dicke = o fo r 



2 



When f2 = 0, eq. ([54]) for T c reduces to eq. (Jl5|) . provided we assume (a Jf )(a) = {a) a) and 
use the relation (a±) 2 = ( a ^)( a )- Thus, both MFT and the coherent state expansion 

lead to the same thermodynamic properties for the Dicke model. 
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t=0, h=0.5 
t=0.12,h=0.5 
t=0, h=0 




Number of Junctions N 



FIG. 1: Coherence order parameter X(N,t), plotted as a function of the number of junctions N, 
for n = and 0.5 at values of the scaled temperature t = and 0.12, as indicated in the legend, 
using J = 0.217, hw = 0.1517, and 5 = 0.1. 



IV. NUMERICAL RESULTS 



We have carried out several illustrative numerical calculations using the models and 
approximations of Sections III Bl and IHIl For the MFT of Section IIIB[ these results are 
obtained by minimizing the mean- field Helmholtz free energy, eq. (1211) . with respect to A 
at each T for fixed g and n. In all our calculations, we have taken hw = 0.15U, J = 
0.2U,g = 0.1, and iV = 110; other parameters are described below. It is convenient to 
introduce a dimensionless temperature t = -fj-- Except for g, these are in the same ratios 
as in recent eperiments of Ref . [25j using the correspondence, (hw, U/8, J) in our notation 
to (hw r = 6.0GHz, Ec ~ 5.0GHz, and Ej^ max ~ 8.0GHz) in the experiment. However, we 
have used a much larger value of g, in order to see the transition to a coherent state at a 
reasonable value of N. 

Figs. [H El [3], HI and [5] show mean-field results for the Hamiltonian of Section [Til including 
all Josephson levels. First, we consider the coherence order parameter \(N, T) assuming g 
independent of N . Fig. [1] shows X(N,T) for n = and 0.5 at t = and 0.12. In all cases, 
there is obviously a threshold number of junctions N c (t) below which A vanishes, and above 
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1 1 ■ L J 1 " u ■ 1 1 

-0.4 -0.2 0.2 0.4 

h 

FIG. 2: Coherence order parameter A(n, t), at several values of the scaled temperature t = 
0,0.12,0.14 and 0.16, as indicated in the legend, using J = 0.2C7, Hiv = 0.15U,g = 0.1, and 
N = 110. 

which A ^ 0. For a sufficiently large N, A — > 1, signaling complete phase locking. Fig. [1] 
shows that both X(t) and iV c (t) decrease with increasing t, and that n = 0.5 leads to a larger 
A at fixed t than does n = 0. Both features are intuitively reasonable, since at n = 0.5, 
the two lowest states of the junction have only a small gap, making it easier to couple the 
junction to the cavity. 

Fig. [2] shows A(n, t) as a function of n, which is related to the voltage across the Josephson 
array, at t — 0, 0.12, 0.14 and 0.16. Since X(n, t) is a periodic function of n with period unity, 
we plot only the range —0.5 < n < 0.5. All the plots of Fig. [2] show that, for any choice of 
the other parameters, A is maximum at n = 0.5. The plots also show that there exist values 
of N such that the array is coherent for some non-zero values of n even if it is incoherent at 
n = 0. Finally, Fig. [2] shows, as expected and as is also shown in Fig. [fl that the effect of 
increasing t at fixed N and g is to suppress A. 

The temperature dependence of A(n, t) is plotted versus t in Fig. [3] for different values 
of n. In all cases, there is a critical temperature t c above which A = 0. Note also that, as 
t — > t c from below, A — ► continuously. This behavior is a hallmark of a continuous phase 
transition. However, since there are only a finite number of junctions, the transition is not 
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0.05 0.1 
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0.15 



FIG. 3: Temperature dependence of coherence order parameter \(t) for n = 0,0.1,0.2,0.3,0.4 and 
0.5, using J = 0.2?7, = 0.15U,g = 0.1, and N = 110. 



a true thermodynamic phase transition. 

We have also calculated the average Cooper pair difference (rij) across the ] th junction, 
within the mean- field approximation. Since all the junctions are assumed identical, (rij) is 
independent of j and may be denoted (n). (n) is related to the voltage drop V across a 
junction by CV/2e = (n), and can be calculated from the relation 

V p -f3E(v=2n+2k,q)( \ 

( n > = ~y p - m u=2n + 2k, q) ■ ( 66 ) 

Z^fc=0,±l,±2,... e 

where (n)^, the expectation value of the operator n in state k, is 
(n)fc = J VC=2fl+2fc(0) y^) ^=2n+2fc(0)# = n+ij y v =2n+2k{v) dyu - 2 '^+ 2k ^ dv. (67) 

The last expression is obtained using the relations n = i-^ = and ip v {4>) = 

e~ m ^~ a ^u u ((j) — a) = e~ 2mv y u (v). In Fig. 0], we show this calculated (n)(n,t) versus n 
for several values of t. Note that, t ^> t c , (n) ~ n. 

Next, we discuss the temperature dependence of the photon probability distribution 
P(n, t) calculated in this mean-field approximation. P(n, t) is defined simply as the probabil- 
ity that the cavity contains exactly n photons at temperature t. Since there is no coherence 
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FIG. 4: Average Cooper pair number difference (n) at t = ksT/U = 0,0.12,0.16,0.2, and 0.24. In 
all cases, J = 0.2*7, fno = 0.15*7, 3 = 0.1, and N = 110 where t c (n = 0) = 0.124 and t c (n = 0.5) = 
0.181. 



for t > t c , P(n) is given simply by the usual Bose distribution with zero chemical potential: 



P(n) 



-/3ftw(n+l/2) 



[eq. m\ 



Using the solutions of H 



E^o e-0Mi+i/2) 2sinh/?^/2 
P(n) = (n\p\n). 



m 

photon ' 



00 

— Y< 

n / j 



p = — 2^ e ~ m e~ ip{x)/h \l) (I I e ip{x)/h 

^photon i = q 



where Z™ hoton is the partition function corresponding to H^ oton . Thus, we obtain, 



PE t / n \ e -ip{x}/h\j\/j\ iD ip{x)/n 



\l)(l\e vp{x)/h \n) = 



(n\e- ip ^l h \l)\ 2 e-f* E K 



1=0 



1=0 



where 



(n\e- ip{x)/h \l) = (n\ 



\x)(x\ 



-ip(x)/h\ 



(68) 



(69) 



(70) 



(71) 



10= / r n {x+(x))Mx)dx. (72) 
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Photon Number 

FIG. 5: Photon number distribution P(n) for the Josephson junction model at various tempera- 
tures: t = 0, 0.12, 0.13, and 0.14. In all cases, J = 0.2*7, huj = 0.15*7, g = 0.1, N = 110, and n = 0. 
For these parameters, t c = 0.131. P(n) represents the probability that there are exactly n photons 
in the cavity mode. 



For the consideration of only I = term, the probability function P(n) corresponds to 

I \2n 

P (n) = ^Ve-H 2 , (73) 



n\ 

where Po(n) has a maximum at n = \a\ 2 . 

Fig. [5] shows this photon distribution at t = fc^X '/U = 0,0.12,0.13, and 0.14 where 
t c = 0.131. From low t up to near t c , P(n,t) is substantial over a wide range of n, but for 
t > t c , the population of the photon state with n = rapidly increases. 

Next, we compare the mean-field results of Section [Til (which includes all junction levels) 
to the results of Section [TTTJ for the Dicke model and the modified Dicke model. We consider 
specifically n = 0.5; at this value of n, the two-level approximation may be best, because the 
two lowest junction levels are separated by the largest gap from the higher levels. In order 
to compare the three models, we plot t c (N) in Fig. El(a), and the average photon number 
(da)(t) in Fig. ©(b). For comparison purposes, we choose the Dicke parameter £ = — 
and the Dicke parameter e = E{2 — 2n,q) — E(—2n,q) = ((7/8) [a(2 — 2n,q) — a(—2n,q)] 



261 ] . Also, we treat Q simply as a parameter determined by best fitting to the results of the 



MFT. We denote the critical number of junctions and the critical temperature of the MFT 
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by N™(t), t™; of the Dicke model, by N^ lcke (t), tf icke ; and of the modified Dicke model, 
by N^ Dtcke (t), t^ IDtcke . The t™'s are obtained numerically; the other t c 's are obtained from 
eq. (I4"5l) for e and from eq. (1641) for t^ Dlcke m i n the mean-field case, when the coherence 
order parameter vanishes we just have the Bose result for the average photon occupation 
number: 

(a!a) m {t > C) = u! • (74) 

e k B T _ I 

On the other hand, in the coherent state, we have 

oo 

(ata> ro (t<t™)=5>P(n). (75) 

n=0 

For the other two models, (a) a) Dicke and (a) a) m Dicke are calculated from the conditions (]4*6l) 
and ([63]) with (ct_l) 2 = (f|) 2 (a t a), respectively 

In Fig. El we plot t c (N) and the average photon number (a"'"a)(t), for the three models. 
To compare the MFT with the modified Dicke model, we have considered four choices for 
Q, corresponding to Qhw/\^\ 2 = 0.2, 0.4, 0.6 and 0.8. All three models show the same 
qualitative behavior, i. e., a transition from coherence to incoherence with decreasing N or 
increasing t. The solutions of the Dicke and modified Dicke models are qualitatively in good 
agreement with that of the MFT; however, the solution of the modified Dicke model with 
Q = 0.8\C,\ 2 /(huj) agrees better with MFT than do any of the other three. 

The behavior of (a^a) (t) differs somewhat among the three models. For the MFT model, 
(a^a)(t) — > as t — > t™ from below, and remains very small, but non-zero, for t > t™. On 
the other hand, in both the Dicke and modified Dicke models, (a^a) reaches exactly zero at 
t = t c and remains zero for t > t c . The most conspicuous qualitative difference between the 
two models occurs at large t, where the Dicke and modified Dicke models give {a^a)(t) = 0, 
while (a j a)(t) in the MFT increases with increasing t according to the Bose distribution. 
This discrepancy probably occurs because the first two models include only two levels per 
junction, while the MFT of Section [Til treats a many-level system. 

V. THERMODYNAMIC LIMIT 

We now make more precise the connection between the Josephson-cavity model and the 
Dicke and modified Dicke models in the thermodynamic limit. We first consider the Dicke 
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FIG. 6: Comparison between the predictions of the mean-field approximation, the Dicke model, and 
the modified Dicke model for the critical temperature [part (a)] and the average photon number, 
(ata)(i) for N = 70 [part (b)], at n = 0.5. For the Josephson junction model, we use the parameters 
J = 0.2£7, hw = 0.1517, and g = 0.1; these lead to iV c m (0) = 64.0 and t™ = 0.0681. The 
corresponding critical numbers for the Dicke model are Nj? lcke {Q) = 37.1 and t® lcke = 0.168. For 
the modified Dicke model, we show plots with VLhw/\i\ 2 = 0.2, 0.4, 0.6 and 0.8. 
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model at T = 0. For N identical two-level systems, the condition for the onset of coherence 
at T = is given by eq. (j44p . With the assumption £ = £/yN, this condition becomes 

-^ = hw. (76) 

To map the Josephson-cavity model onto the Dicke model, we assume that the Josephson 
coupling parameter g oc 1/ \/N. As mentioned earlier, this assumption seems reasonable in 
the thermodynamic limit, since according to eq. ([6]), g oc 1/W for fixed cavity shape. 
Writing g = g/y/N, we may express the Josephson coherence condition (|23|) as 

|a'(0)|rf J = /^>. (77) 

We now show that this condition reduces to eq. (I76j) in the limit J < [/ '. The eigenvalue 
derivative |o'(0) | can be obtained approximately in this limit by differentiating the right-hand 
side of eq. (fT8l) with respect to q. Substituting back into eq. (J77I) . we obtain 

, m? = (78) 

v/f/ 2 (l-4n 2 ) 2 + 4J 2 

We can also compute the splitting between the ground and first excited states. If < n < 
1/2, and J <^ U, it is easily shown that the splitting AE between the ground and first 



excited states of the Josephson junction is equal AE ~ \ [U 2 (l — 2n) 2 + AJ 2 ] 1 ^ 2 27j]. Using 



the approximation, 1 — 2n ~ 1 — 4n 2 , we can rewrite the coherence condition ( 1751) as 

-y- = ho. (79) 

This condition is identical to eq. f[TB"j) . with the identification e <-> AE 1 , £ <-> The 
parameter f2 has primarily a quantitative effect on the coherence transition in this model. 
As discussed earlier, in order for the modified Dicke model to be well-behaved in the ther- 
modynamic limit, Q must vary as 1/N. Therefore, we write Q = 2Q/(N — 1); we use — 1 
rather than since each two-level system interacts with — 1 others. Substituting this 
relation into eq. (1601) . and again using £ = ^/x/TV, we find that eq. ( 1761) is replaced by 

4 «" 2 ; to "> = (so 

If the left-hand side is larger than hui, the system is coherent at T = in the thermodynamic 
limit; otherwise, it is not. Thus, a positive Q actually inhibits coherence at T = (not 
unexpectedly, since a positive fl represents a repulsive interaction). 
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Similarly, one can recalculate kBT^ IDtcke [eq. (1M|) using the above iV-dependence of £ and 
Q, with the result 

k B T c MDlcke = i r . (81) 



2 tanh 



1 / hjje 



4^2 ^ hujCl 



Thus, for given values of u, e, and £, the coherence transition temperature is reduced by a 
finite Cl, showing that this form of direct interaction between junctions inhibits coherence 
in the modified Dicke model. 



VI. DISCUSSION 



In this paper, we have calculated the equilibrium properties of an array of identical 
Josephson junctions coupled to a single-mode electromagnetic cavity at temperature T, by 
generalizing a T = MFT [9]. Within the MFT, this system shows a continuous transition 
between coherence and incoherence at a critical temperature T c , provided that the number 
of junctions N > N c . We have also compared our mean-field results to the solutions of the 
Dicke and modified Dicke models. When the parameters of the Dicke model are adjusted to 
match those of the Josephson- cavity system, the two approaches agree qualitatively. 

Next, we briefly discuss the expected accuracy of our mean-field approach, used in Sec- 
tions [Til and Hill The MFT appears reasonable, because for both models, all the junctions 
(or all the two-level systems) interact with the same harmonic mode and hence, in effect, 
with all the other junctions or two-level systems. Since each junction or two-level system 
effectively has many "neighbors", there are only small fluctuations in the environment of 
each about its mean, provided that N is sufficiently large. Thus, the mean-field approach 
should work well at large N. In support of this picture, we have shown that, when the 
mean-field approach is applied to the Dicke model, it produces the exact result (obtained 
from a coherent state expansion). 

Besides the mean-field approximation in Section [Til we have expanded the Josephson 
coupling in powers of the interaction parameter g(a + a)). The value of this quantity can 
be estimated as follows. g(a + a)) ~ g(a + a^) = 2Jg 2 NX/ (huj) = 2Jg 2 X/ {ftw) if we assume 
g = g/\fN. Since A -C 1 is small near T c , this approximation is accurate in this regime, but 
may break down deep in the coherent regime. Therefore, a more accurate approach than 
that used in Section [IT] may be desirable in order to treat the entire regime < t < t c . 
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We briefly comment on the nature of the coherence transition emerging from our mean- 
field approach. This approach produces a continuous transition, i.e., the coherence order 
parameter A varies continuously with t. By contrast, another recent calculation finds a 
first-order transition, in which there is a discontinuous jump in the order parameter at the 
superradiant transition. Their Hamiltonian also has the form of a generalized Dicke model, 
but slightly different from ours: 

N f A 1 

H lee = ata + £ — =(a + <J){a{ + at) + -a{ - Ja^ 1 } , (82) 

3=1 

where a 3 +1 o~j, and a{ are the usual Pauli spin operators for the ] th two-level system. Our 
mean-field treatment of our own generalized Dicke model does not give a first-order transi- 
tion. We speculate that the difference is due to a real distinction between the two models: 
Hi ee has only nearest-neighbor interactions between spins, in addition to the usual Dicke-type 
model, whereas our modified Dicke Hamiltonian has an additional term which is long-range. 
Perhaps the long-range nature of this additional term helps to maintain the continuous 
nature of the coherence transition, as well as the accuracy of MFT. 

Finally, we discuss how our model could be generalized in order to make it more a more 
realistic basis for treating Josephson arrays in a cavity. For real systems, there are other 
factors affecting Josephson junctions besides those included here. For example, there are 
effects due to dissipation, either due to the finite Q of the cavity, or a finite dissipation within 
individual junctions. Both effects can be treated by considering the Josephson junctions as 
coupled to appropriate baths of harmonic oscillators 29, 3o|. When these dissipative degrees 



of freedom are properly included, the nature of the coherence transition may be changed. We 
plan to include some of these effects, as well as the effects of disorder, in a future publication. 
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